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The Venus Balloon Project Telemetry Processing 

J. M. Urech, A. Chamarro, J. L. Morales, and M. A. Urech 
Madrid Deep Space Communications Complex 


j The peculiarities of the Venus Balloon telemetry system required the development of 
j a new methodology for the telemetry processing ; since the capabilities of the DSN telem- 

j etry system do not include burst processing of short frames with two different bit rates 

j and first bit acquisition. The MDSCC engineering staff was tasked to produce a software 
package for the non-real time detection , demodulation, and decoding of the telemetry 
streams obtained from an open loop recording utilizing the DSN spectrum processing 
subsystem-radio science (DSP-RS). ZSft urticle^contains m general description of the 
resulting software package ( DMO-5539-SP ) and its adaptability to the real mission's 
| . variations f ^ CO\ o'* i .< "( * 


I. Introduction 

The Venus Balloon telemetry consisted of 330 seconds of 
1.67-GHz emissions every 30 minutes, one hour, or two hours, 
depending on the mission encoder phase. The 330-second 
emission sequence began with 30 seconds of carrier only (TM1 
mode), followed by 270 seconds of telemetry data (TM2 
mode), and ended with another 30-second period of carrier 
only. 

During the first 210 seconds, the telemetry data rate was 
4 bits/sec, and the last 60 seconds, the rate was 1 bit/sec (see 
Fig. 1). This telemetry data was convolutionally coded ( K = 6, 
V= 1/2), and thereafter Manchester encoded with a subcarrier 
derived from the carrier (1667.9 MHz) divided by 6,553,600, 
resulting in an approximate subcarrier of 254.5 Hz (see Fig. 2). 

The DSN spectrum processing-radio science (DSP-RS) 
recorded 20K samples/sec sampled spectra on magnetic tapes 
(level 0). These tapes were processed by DMO-5539-SP (Venus 
Balloon telemetry software), thus obtaining the decoded data 
(level 1) tapes. 


II. Telemetry Software Description 

The telemetry software package ran in a DSN spectrum pro- 
cessing (DSP) environment with three additional peripherals: a 
local terminal, a hard copy device, and a color display. These 
peripherals enabled a stand-alone operation isolated from 
other station activities. 

A. Design Criteria 

Design parameters were given as (a) signal strength: nominal 
18 dB (1 Hz bandwidth), with uncertainty as low as 13 dB; 
and (b) 1.67 GHz carrier stability: 


The 1.67-GHz carrier uncertainty was estimated in ±2 kHz, 
Doppler rate 0.05 Hz/sec. Based on these parameters, the 
design was concentrated on the detection of weak signals 
highly stable in frequency and affected by moderate Doppler 
variations. 


199 



Due to the carrier uncertainty, the bandwidth of the radio 
science open-loop receiver (MMR) analog filter selected was 
4091 Hz for frames A-2 to A-6, B-2 to B-46, and A-l and B-l 
(the first frames of every Balloon); when uncertainty was higher, 
the filter selected was a 8182-Hz bandpass filter. The first 
Balloon was designated as “A” and the second as “B.” Teleme- 
try frames for each Balloon were numbered 1 through 46. 

With these constraints in mind, the software was designed 
with sophisticated algorithms to confirm the presence of the 
1.67-GHz carrier under noisy conditions and proved to work 
satisfactorily with 6 dB less than the nominal case during the 
compatibility tests with the spacecraft prototype test simula- 
tor (PTS). During the real mission, the received signal strength 
was close to nominal, but the frequency stability and Doppler 
variations were worse than anticipated by several orders of 
magnitude. The software was highly modular, and its adapt- 
ability to the real situations was relatively simple, although the 
processing of several “pathological” frames required some 
interactive operation. 

B. Preprocessing 

The DSP-RS recorded tapes were selectively transferred to a 
large capacity disk (LCD) from which the spectrum data was 
further processed. During this transfer, only 350 seconds of 
recorded data were copied from the magnetic tape onto the 
LCDs containing the 330-second emission of the Balloon. 

This transfer also performed a 1/2 decimation, reducing the 
20K samples/sec to 10K samples/sec. The determination of the 
MMR filter shape was also possible at this stage. The resulting 
volume of data to be processed was 3.5 millions of samples 
(8 bits wide) per Balloon emission. 

C. Telemetry Processing 

A typical process consisted of the following major steps: 

(1) Carrier confirmation 

(2) Carrier estimation 

(3) Doppler estimation 

(4) Carrier demodulation 

(5) Subcarrier demodulation 

(6) Symbol and frame synchronization 

(7) Viterbi decoding 

1. Carrier Confirmation. The purpose of this step was to 
determine whether a Balloon telemetry frame existed in the 
recorded data or not. To this end, the carrier confirmation 
module used 30 seconds of recorded data starting from a 


hypothesized carrier-only period, which is to occur in the 
beginning of each telemetry frame (see Fig. 1). 

The process was initiated computing 4096-point FFTs 
every 0.5 sec of recorded data, thus obtaining 60 power 
spectra. Then the ten maxima method is applied, forming a 
histogram from the 60 sets of 10 maxima of each spectrum. 
In this histogram, the confirmation module may find a zone 
where n consecutive bins contain a cumulative occurrence of 
maxima. After an ascending time search of the 60 power 
spectra, a maxima belonging to that zone satisfying a thresh- 
old value v is a candidate for the searched carrier bin. If m 
consecutive power spectra contain maxima in the same zone, 
the presence of carrier is confirmed; if not, the first power 
spectrum of the set of 60 is discarded, the hypothesized time 
start is delayed by 0.5 sec, a new FFT is computed for the 
next 0.5 sec of recorded data, and a new set of 60 power 
spectra is analyzed with the maxima method. This process 
is repeated until either the carrier is confirmed or a recording 
time limit is reached. Figures 3 through 6 illustrate a nominal 
case and weak signal plots. (Typical values of parameters n, v, 
andm were 10, 15 and 2, respectively.) 

Once the carrier start time has been determined, the maxi- 
mum of each of the 60 FFTs is passed to the following module. 

2. Carrier Estimation. In this step, the software performs 
the initial estimation of the carrier frequency and rate using 
the maxima of the previous step, whenever they correspond to 
the zone of bins determined there. 

Bin numbers are converted to Hertz, thus obtaining the 
translated 1 .66-Hz carrier frequency referenced to the origin 
of the MMR filter used to condition the DSP-RS discretization 
process. These frequencies are fitted to a straight line using the 
least squares method. The resulting regression line provides 
then the initial estimations of carrier frequency and rate dur- 
ing the first 30 seconds of the Balloon telemetry emission. If 
the error of a particular point (i.e., its distance to the regres- 
sion line) is greater than a delta / parameter (typically equal 
to 5 Hz), then the point is eliminated, and the fitting is 
repeated. 

3. Doppler Estimation. After the initial estimation of 
carrier frequency and rate obtained in the first 30 seconds of 
the telemetry transmission by the previous module, the Dop- 
pler estimation module worked over the entire 330-second 
period of the Balloon emission. 

In this case, the software algorithm was based on a least 
squares non-linear regression process with the following 
model: 
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M.it) = A.sm(y.+f.t+^y.t 2 ) 

The model was applied over spans of recorded data with 
lengths of 10, 4, and 2 seconds, respectively, chosen depend- 
ing on the observed carrier stability and Doppler rate. In 
addition to a complete run over the 330 second period, the 
software was able to be directed to repeat the process on 
individual spans or groups of variable numbers of spans, with 
different span lengths or initial parameters. 

Each span overlapped the precedent one by 50 percent. 
The stepping from one span to the next one was based on the 
SNR of the fitted model; if it violated a pre-established thresh- 
old, the model was iterated with a spiral search pattern until 
the SNR threshold was satisfied. The whole process is detailed 
in Appendix A, and illustrated in Figs. 7 through 10. 

4. Carrier Demodulation. The carrier estimated by the 
previous module is demodulated from the recorded data 
simples, S., using the following method to maintain phase 
continuity. 

S. = S. [(1 - X/) cos co t (f.) + Xz cos co (f.)] 

i i LV ' n- 1 v i' n K i /J 

where 

r.- t 

X . = 1 _ ” — , n - 1 < i < n 

l n n—1 

co n (0 = carrier estimate in span ri 

The DC component is removed before and after the demodula- 
tion in groups of 5000 samples. After this process, the demodu- 
lation samples are reduced from 10K samples/ sec to approxi- 
mately 2500/samples/sec, using a running averages method. 
The resulting bandwidth is then approximately 1250 Hz. 

5. Subcarrier Estimation. The first step in this module is to 
establish the coherency between the spacecraft timing (in sub- 
carrier cycles) and the ground station timing, from which the 
data has been synchronized during the sampling and recording 
by the DSP-RS. Computing the received carrier from the open- 
loop receiver local oscillator frequency and the translated 
signal in the MMR filter, we then obtain the actual spacecraft 
subcarrier frequency by dividing the detected carrier fre- 
quency by 6,553,600. 

After this frequency is known, the sampled data is arranged 
so that each subcarrier cycle contains 10 samples, as follows: 


s s Ci = j-f B(t)dt 

sJt i 

where 

S . = subcarrier sample 

B(t) = modulated subcarrier 

t _ sample rate 

s 10 X subcarrier frequency 

From this point on, every pseudosymbol in the last 60 seconds 
of the telemetry frame contains 640 subcarrier samples, and 
in the first 210 seconds contains 160 subcarrier samples. 

The software now focuses on the last 60 seconds of the 
telemetry frame, and to estimate the subcarrier phase, uses 
five subcarrier models with a relative delay of one tenth of 
a pseudosymbol (i.e., 64 subcarrier samples) and computes 
the correlation function of the five models in groups of 64 
subcarrier samples over the 60 second period of 1 bit/sec (or 
4 pseudosymbols/sec) (see Fig. 11). The model of maximum 
correlation determines the subcarrier phase delay in samples 
and a first estimation of the pseudosymbol start in tenths of 
pseudosymbols. 

6. Subcarrier Demodulation. The demodulation process is 
accomplished by taking groups of 10 subcarrier samples (one 
subcarrier cycle) and inverting the last five, using the subcar- 
rier phase information provided by the previous module. At 
the same time, a bandwidth reduction is performed averaging 
the 10-sample groups. The resulting bandwidth is 128 Hz. 

7. Symbol Synchronization. This module continues to use 
the 60 seconds of 1 bit/sec (4 pseudosymbols/sec) to deter- 
mine the symbol synchronization. The method is based on 
averaging two half-pseudosymbols (i.e., two groups of 32 
samples) centered on the estimated crossing from a positive 
pseudosymbol to the adjacent negative pseudosymbol. 

The difference between the two averages is proportional to 
the crossing estimation error, which is iteratively corrected 
until the crossing error is less than or equal to five samples. 
This process also keeps track of sign transitions to determine 
the correct pseudosymbol pairing (the two pseudosymbols of 
the same symbol must have opposite signs.) 

The final estimation of the symbol start delay is obtained 
after computing the pseudosymbol SNRs with a variation of 
up to 24 samples from the last crossing estimation. 
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When this symbol synchronization process is completed, 
the frame is normalized to 16 samples per second (this timing 
is referred to the spacecraft timing using the determination of 
the subcarrier coherency in subsection 5). The frame has been 
now reduced to 5280 samples. 

8. Frame Synchronization. The Balloon telemetry frame 
consists of 900 bits with the following distribution: 

48 bits (synch pattern) transmitted at 4 bits/sec 

792 bits from 22 measurements of 6 instruments (6 bits 

each), transmitted at 4 bits/sec 

60 bits of 10 data channels (6 bits each), transmitted at 

1 bit/sec 

The 96 symbols of the frame synch information are contained 
in 384 samples of the reduced set. This software module looks 
for the synch pattern correlating the sample set starting at 
500 different sample positions around the expected point. The 
maximum correlation indicates the position of the synch word 
in the sample set (see Figs. 12 and 13). 

From this position, the 1680 symbols corresponding to the 
4 bits/sec are computed by removing the Manchester encoding. 
The 120 symbols corresponding to the 1 bit/sec are obtained 
in a similar fashion from the remaining samples, averaging four 
samples per pseudosymbol. 

9. Viterbi Decoding. This module performs the convolu- 
tional decoding of the 1800-symbol frame using a standard 
algorithm. 


After completing the decoding, the resulting bits are con- 
volutionally encoded to obtain a symbol reference frame, 
which in turn is used to compute the symbol error distribu- 
tion. This is presented together with the decoded frame as 
the end of the overall process. 

The analysis of the synch word pattern permits the detec- 
tion of a data bar condition. Final results are given in data 
that are true and decommutated (see Fig. 14). 


III. Adaptability to the Real Conditions 

During the mission, the spacecraft received signal pre- 
sented the following characteristics: 

(1) Average carrier SNR in a 1-Hz bandwidth approxi- 
mately 20 dB (expected value 18 dB) 

(2) Average symbol SNR approximately 7.5 dB (expected 
6.0 dB) 

(3) Frequency variations expected 0.05 Hz/sec 

(4) Frequency variations received up to 3.0 Hz/sec 

Frames processed with the antenna near horizon and toward 
the end of spacecraft battery’s lifetime have lower SNR. 

The software was able to cope with the above frequency 
variations after creating new configurations with existing or 
modified modules and optimizing Doppler estimation param- 
eters. These software modifications resulted in an average 
degradation of 0.5 dB. 
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Fig. 1. Telemetry mode sequence 
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Fig. 2. Telemetry modulation scheme 
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Fig. 3. Carrier detection FFT, nominal signal 


Fig. 5. Carrier detection FFT, weak signal 
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Fig. 4. Carrier detection histogram, nominal signal 
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Fig. 6. Carrier detection histogram, weak signal 







Fig. 8. Carrier frequency vs time (Balloon A2) 


TIME 162 

‘ 4”23 ”45.80 

80 

PHASE 

J IT. b H~ 

EKEU. 

I 4 b 4 . 1 b 

KATE -2.6821 

amp. 

0 . 0 88 

SN ft 

1 / .00 

TIME" 167” 

' 4 Zj 'ifj.aO SPATE 

6 1 

Pri Ao£‘ 

29T.11 

r*KE0. 

"14 5 0.50 

HATE - U . 4 8 b S 

AMP. 

0.092 

SNK 

1 / . b 2 

TIME 162" 

4* 2 8 *r # 8‘0 “SPAM ” 

8 2 

PHASE ' 

HI .29 

FKEQ.” 

1448.77 

KATE 2.110b 

AMP. 

0.100 

sriK 

1 / .99 

TIME i 62 

4 2 i 4 3. HO "SPA iT” 

6 S 

PHASE 

2b. /9 

FKEvJ • 

1450.40 

KATE 8.2802 

AMP. 

0 • i 2 8 

SN ft 

2i.il 

TIME 162 

4 21 V/.HG'SPA.r 

8 4 

PHASE ' 

vrr/- 

Fkeu. 

1454.28 

KATE 2 . U 2 1 9 

AMP. 

0.09 / 

SN K 

14.68 

TIME' r 62 

- 50.<ro “SPAT J — 

' 8 5 

Pn A 3E 

- 2rvitf 

EKE U . 

14b 7.71 

KATE - 0 . 7 9 o 4 

AMP. 

0.101 

SNft 

14. / 2 

TIME 162 

V”2“8 _ bl . 80 "3P7UT” 

"Ob 

PHASE 

bj ,9b 

Ekeu . 

1 4 b 8 . t u 

KATE -8.1715 

AMP. 

U . 082 

SN ft 

16.22 

T r M E "" 1 67” 

*’4“2i"52.30"“SPA7;“ 

87 

PnAoE “ 

216. b6 

FftEU.' 

1455. 2 1" 

KATE -8 .7863 

AMP. 

0.092 

b N ft 

Id. 84 

TIME' r62 

4 2 J b J.3G 5 PA M 

88 

PHASE - 

“1 6 • b 8 " 

EKEu ; _ 

1 4 50 . b 2 

KATE -2.0800 

AMP . 

0.089 

SNH 

l/.il 

TIME 1 6Z ” 

4 71‘ 54;30 oPAM-” 

“6 9 

PHASE 

282.68" 

FKEO. - 

1447.16 

KATE " 0.246b 

AMP. 

0.0/6 

SN ft 

1 6 • 1 4 


Fig. 9. Partial listing of Doppler estimation (Balloon A4) 
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Fig. 12. Frame synch word detection 
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Fig. 13. Decommutated frame 




Appendix A 

Doppler Estimation Span Processing 


The received signal (PCM/PSK/PM) can be represented 
mathematically by (Refs. A-l , A-2): 

S(t) - y/2P (sin 00(f) cos 0 + cos co(f) sin 0 B(t)) 

This signal is mixed with a local estimate: 

S m (t) = cos w(f) 

and the result, after a low pass filter provides S f samples that 
are a function of the difference oo - co. The process is as 
follows. 

Mixing ( see Fig. A-l): 


o> D (t) = <o(0- <£(r) = \p +ft +jyt 2 

The problem now is to minimize 

E {S r M) 2 

Since S. is equal to the reduced sample set 

M. = A sin (u D (t.)) = A sin (i// Aft. + 7 t?} 

that is equivalent to solving the non-linear system: 


S(t) • S m ( t ) = y/2P (sin 00(f) cos 0 + cos co(f) sin 0 B(t)) 

X cose o(t) 

= \/27> (sin (c o(f) - ou(f)) cos 0 + ^ 

X (sin (co(f) + co(f)) cos 0J 
+ \f2P ((cos (co(f)) - L o(f) + (cos (00(f) 

+ eo(f))) sin 0 £(f)J 

Filtering: 


V' 7*™/ 

- 22 .( 5 .-^.)^= 0 

where 

1 < k < 4 



Expanding this expression, we obtain the system in explicit 
form: 


f+Af/2 


w 

Jt- At/2 


This filter cancels terms with 
B{t), resulting in: 


5(r) • SJt)dt 

(co(t) + <o(f) and terms with 


5. = 
2 


y/2? 


sin [co(f ) - 00(f)] cos 6 


and calling: 


A 


s/2P 

~ T 


cos 0 


i 

i 

EiV ' 1 

i 


sin 0J d (r, )) sin u d {t.) = 0 


sinoo d (f.)>4 cosoo rf (f.) = 0 

sin u> d (t.)) At cos cj d (t.) = 0 

At 2 

sin co d (t.))— cos w d (t.) = 0 
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After sample trigonometric operations, the system is 
expressed with the following non-linear functions: sin co d , 
cos u> d , A sin 2co d , A cos 2co d . This non-linear system is 
linearized to approximate the above functions by their Taylor 
series expansions. 

After establishing a linear system of four equations with 
four unknowns (i//, /, 7, A ), the solution is obtained by itera- 
tion from a starting point (^ 0 ,/ 0 , 7 0 M 0 )- Iteration incre- 
ments for / and 7 are obtained in a spiral way (see Fig. A-2), 
where negative values of / are avoided. 

Limits of the spiral way are: 


The iteration process is started with an initial vector 

^0 ~ (^0 ’^0 ’ ) 

and an incremental vector: 

4*0 = (^ 0 > A ^ A vM>) 

in the usual way: 

= *o +A *o^ A *i 


where 


X N = X N - 1 +A ^V-1 


AX 


TV 


The process finishes when: 


/ = span length in seconds 


lAX^KlO" 5 , 1<*<4 


BW = low pass filter bandwidth 

Typical values are: f b = 0.8 Hz,^ im =0.8 Hz,7 lim =0.05 
Hz/sec. The initial values for the other two unknowns are: 

i 

(initial amplitude) 

\h = K \b 

with K such that it maximizes: 

£s,*in(K* 10 */ 0 ',+ ‘ V?) 

i 

for 1 < K < 18 and^ 20 = phase corresponding to the first step 
of 20 degrees. 


If this limit is not reached, the sequence is aborted when N 
is greater than the other limit (typically 60). Then, new initial 
values are obtained from the spiral search and the process is 
resumed. 

If the convergence to 10~ 5 is obtained, a carrier SNR is 
computed and tested against a threshold value. If the test 
passes, this vector (v^,/, 7, A ) is considered a valid solution. If 
the test fails, the process is repeated with new initial values 
from the spiral search until a solution is found, or the spiral 
way is exhausted. When the latter happens, the solution with 
higher SNR is adopted. 

Multiplicity of Solutions : If (1 p, /, 7, A) is a solution, 
(-i//, -7, -A) and (\p +N , -/, -7,^4) are possible solu- 

tions that must be discarded. The set (-1//, -/, — 7, -A) is easily 
detected by checking the sine of the amplitude (X). The set 
(\p +N , -/, -7, A) is avoided by mixing with a frequency lower 
than the one estimated (typically 0.8 Hz lower), so that after 
mixing the difference will always be positive. 
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Fig. A-1. Doppler estimation span processing 
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